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Abstract 



Aims. On the basis the theoretical model proposed by Bekenstein for a's variation, we analyze the equations that 
describe the energy exchange between matter and both the electromagnetic and the scalar fields. 

Methods. We determine how the energy fiow of the material is modified by the presence of a scalar field. We estimate 
the total magnetic energy of matter from the "sum rules techniques" . We compare the results with data obtained from 
the thermal evolution of the Earth and other planets. 

Results. We obtain stringent upper limits to the variations in a that are comparable with those obtained from atomic 
clock frequency variations. 

Conclusions. Our constraints imply that the fundamental length scale of Bekenstein's theory "^s" cannot be larger than 
Planck's length "^p". 

Key words. Bekenstein's model-planetary heat flux 



1. Introduction 

The time variation in the fine structure constant has been 
studie d several times since first being proposed bv lGamowl 
()1967() . Observational upper bounds on its time variation as 
well as several theoretical frameworks that consider a as a 
dynamical fi eld have been published (a n exhaustive list can 
be found in ()Landaull2002l : IUzanll2003D and references there 
in). Although still disputed, the clai m that a was smaller in 
the past is an exciting pe rspective, (iMurphv et a l. 2003). 

Beckenstein's theory ()Bekensteinlll982l) . which is based 
on a number of minimal hypothesis of highly accepted phys- 
ical principles, is in a sense representative of many low en- 
ergy theories inspired by grand unification schemes. In this 
work, wc derive equations that govern the energy exchange 
between matter, the scalar field, and the electromagnetic 
field. Although we do not analyze the precise mechanism of 
energy release, we assume that the work done by the scalar 
field is radiated away in an efiicient way, as for the roto- 
chcmical hea ting of neutron st a rs due to the spin down 
of th e star (jReiseneggeil Il995t [Fernandez fc Reiseneggerl 
I2005D . 

In section [21 we briefiy review Beckenstein's theory, as 
well as the cosmological time evolution of a that it predicts. 
In section[3l we derive a generalized version of the Poynting 



* fellow of CONICET 

** member of the Carrera del Investigador Cienti'fico y 
Tecnologico, CONICET 



theorem for the electromagnetic field, and from the conser- 
vation of the total energy-momentum tensor we find how 
the energy flow of matter is modified by the scalar field. In 
section |4l we discuss the magnetic energy of matter using 
a simple nuclear model. In section [SJ we study the thermal 
history of the Earth in the presence of Bekenstein's scalar 
field. We also describe in section |6] the results we obtained 
for the outer planets. Finally in section [7] we summarize our 
conclusions. 



2. Time variation of a in Bekenstein's formalism 

We briefiy review Bekenstein's formalism and its predic- 
tion for the cosmological time variation of a. Although we 
consider galactic as well as terrestrial phenomena, we can 
nevertheless confidentl y assume that they tra ck the cosmo- 
log ical evolution of a , (jShaw fc Barrowl[2006l) . 

iBekensteh] (|1982D proposes to modify Maxwell's theory 
by introducing a field e that dynamically describes the vari- 
ation of a. The foundational hy pothesis are the following 
(|Bekensteinlll982t lLandaull2002() : 

1. The theory must reduce to Maxwell's when a = Cte. 

2. The changes in a are dynamical (i.e. generated by a 
dynamical field e). 

3. The dynamics of the electromagnetic field, as well as 
the e field can be obtained from a variational principle. 

4. The theory must be locally gauge invariant. 
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5. The theory must preserve causahty. 

6. The action must be time reversal invariant. 

7. Planck's scale is the smallest length available in the 
theory. 

8. Einstein's equations describe gravitation. 

String theories and the like in which there are other 
fundamental length scales, force us to set aside condition 
[71 These hypothesis uniquely lead to the action 



where 

Se = 



1 

16^ 

he 

21^ 



S = S, 



f 5, + 5, 



Sc 



(1) 

(2) 
(3) 



Sm and Sg are the matter and gravitational field actions, 
respectively, is the so-called Bekenstein's fundamental 
length, and the metric here is (—1,1,1,1). 

The main difference between Maxwell's and 
Bekenstein's theories is the connection between the 
vector potential and the electromagnetic field 



1 



(4) 



and the second kind of local gauge invariance implies that 



eA,. 



X., 



(5) 
(6) 



as the gauge transformation and covariant derivative of 
the theory, respectively. The last equation defines the local 
value of the elementary electric charge (coupling constant) 



e{r,t) = eoe(r,t). 



that is 



(7) 



(8) 



In what follows, we neglect the small spatial variations in a 
and focus on the cosmological variation, as we are interested 
in any secular energy injection of the scalar field on a planet 
such as the Earth. In our approximate analysis it is also 
enough to work in a flat space-time. 

The field equations for the electromagnetic field and for 
e are 



e 



(2 

□ lne= ^ 

he 



de 47r ^ ^ 

da 



he V de 



(9a) 



(9b) 



where = Y.{eo/ c-1)u^'{~gy^l'^5^[x' - x\t% uf is an 
"standard estimate" of the 4th velocity of each particle ac- 
cording to the mo del, and a is the energy density of matter, 
iBekensteinI (|1982[ ). □ is the covariant flat d'Alambertian 



□0 



(10) 



A note regarding the matter Lagrangian is in order: in 
(jBekensteinl I1982L l200i) . Bekenstein represents matter as 
an ensemble of classical particles. However, wherever quan- 
tum phenomena become important, as in white dwarfs or 
condensed matter physics, this is not a realistic description. 
It is also an inaccurate description at high energies (or on 
small length scales) because fermions have a "natural length 
scale" , the particle Compton wave-length Ac = ft/mc, 
which makes quite unrealistic any classical mod el at higher 
energi es. In particular, several conclusions of IBekensteinI 
(|2002D have to be reconsidered. 

In IBekensteinI (|1982l ). it is shown that the cosmological 
equation of motion for e is 



dt\e 



he 



da 1 . n 
'Te - 



(11) 



In the non-relativistic regime, ^ and cr oc e^, hence 



dt\ e 



£2 

nc 



(12) 



where pm is the total rest-mass density of electromagnet- 
ically interacting matter and Cc is a parameter describing 
its "electromagnetic content" , which is essentially the frac- 
tional contribution of the electromagn etic energy to the r est 
mass. A first estimation according to IBekensteinI (|2002t) is 



Cc ^ 1.2 X 10 



-3 



(13) 



Following the standard cosmological model, we assume that 
dark matter is to be electromagnetically neutral. 

Given that pm oc a~^, we can integrate Ea. (jl2[) obtain- 
ing 

e 



<c 



which can be written, using the usual cosmological nota- 
tion, as follows 



Stt 



ait) 



(t-tc). (14) 



The primordial nucleosynthesis standard model tells us that 
the integration constant tc must be very small in order 
not to spoil the agreement between theory and observation. 
Using WMAP values, we obtain the prediction for {a/a)o 



Any measurement with a precision such as a{a./ Hqq) ^ 
10"^ is difficult to achieve, so the comparison between the- 
ory and experiment is a difficult task. 

The same arguments can be applied to m any the- 
ories with varying a, such as Kaluza-Klein (jLandaul 
l2002t) or string spired theori es such as Damour-Polyakov's 
(|Damour fc Polvakov|[T994a] [g). 



3. Energy transfer in Bekenstein's formalism 

We study how energy is injected and then released in vary- 
ing a theories, in order to look for observable consequences 
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in the emission of astrophysical as well as geophysical sys- 
tems. The energy momentum tensor in Bekenstein's theory 
is defined as usual to be 



T^"' ^ 2c- 



6S 



(16) 



Using c = 1, the electromagnetic contribution then has the 
same form as in Maxwell's theory 



rjicm 1 



TP IP ^ _ TP TP"^^ 



(17) 



the difference lying in the connection between the vector 
potential and the field given in Eq.Q. 

The energy-momentum tensor of the scalar field e is 



(18) 



^ II V 2^ 

In what follows, we use the redefined field 

V' = lne. (19) 

As we consider local phenomena, we can work in a lo- 
cally inertial coordinate system. We denote the "field part 
of the energy-momentum tensor" as the scalar plus electro- 
magnetic energy momentum tensor 



~ cm + c 



In terms of ip and replacing g'^'^ with rj'^'^ , we obtain 



(20) 



(21) 



The divergence of Tf is 



The equations of motion Eqs.® are 

which can be used in Eq. ip^ obtaining 
- l^^F^^F^p^,] 

This expression can be simplified using the homoge- 
neous Maxwell equations 



(24) 



(25) 



which cancels out the first bracket. The first and last term 
in the second bracket also cancel out, thus we obtain for 
Ea.(P^ the expression 

Tr,, = - (^^rF'^o. 

f da 1 A (26) 

We add to both sides of the equation the divergence 
of the energy momentum tensor of matter ^ in order 
to find the energy transfer (according to the hypothesis 8 
in Sect. [2l we assume that Einstein's equations hold un- 
modified for the gravitational field and hence that the total 
energy momentum tensor is conserved) 



~ ^ra . V ^ J a 



(27) 



This equation explicitly shows the energy transfer from the 
field e to matter 



m M ~ 



which is the source of any observable effect. From 

e^^ 1 a^i, 
~7 ~ 2~' 



pap 



(28) 



(29) 



we find that the "machian" contribution to energy transfer 
is given by 



rpjiv (machian) 1 '^,1' I . rpf^u p P"/* 

~ 2 a ™^ IGtt 

(30) 

Using Bekenstein's notation, that is, if the time-space com- 
ponents of e^F^" are identified with E while space-space 
components are identified with B, the contribution then 
takes the form 



.do 



(22) % 

(23a) 
(23b) 



[machian) 



dijj 

Stt 
da 



(B2 - 



0%! 47r 



(31) 

where S = ^j^. Then, the component of Eq.lgHl) reads 
= j.E - e-^^'^ - e-2^VV..S -f (32) 



47r 



di}} ' 



An implicit assumption of our previous analysis and al- 
gebra is the generalized Poynting theorem. In its standard 
version, it involves only electromagnetic terms, while in our 
case it also involves the interaction between the electromag- 
netic and scalar fields given by 



1 



(33) 
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Using again Eq.d^H), 



em u 



16n 



167r 



em.u 



P aPO'P 

167r 



Then, 



-2^S.VV'-^^(B2-E2), 



Op 



9Ue 



dt 
= -Ej 



V.e- 



Stt 



47r 



47r 



(34) 



(35) 



(36) 



(37) 



where T^m ,o = {duem)/dt, the electromagnetic energy is 
= e-'^'{E' + B2)/(8^), Te™°\, = S/.e-^^{^) = 
V.e~^'''S, and S is the Poynting vector. We note that this 
result is independent of the details of the gravitational and 
matter Lagrangians, as well as their interacting terms with 
the electomagnetic field. In particular, it holds indepen- 
dently of the details of the interaction of matter with the 
scalar field. We recall that the usual interpretation of the 
first term in the right hand side of Eq. ([57|) is the work done 
by the electromagnetic field on matter. In the same fash- 
ion, we may interpret the second and last term as the work 
done by the electromagnetic field on the scalar field. An 
analog phenomenon would be that given by the work done 
by an increasing Newton constant G on a planet augm ents 
its pressure and thus compresses it (jJofre et al.ll2006() . 

We estimate the electrostatic contribution to the matter 
energy. In a non-relativistic system such as a light atom or 
nuclei, the electromagnetic energy is given by the electro- 
static field that satisfies the equation 



(38) 



where is the reference charge density. In the limit where 
a varies only cosmologically we have 



(39) 



whose solution is 



E^e^'f'Eo, (40) 

where Eq is the electrostatic reference field defined for e'^ = 
1. We can evaluate the electromagnetic energy density to 
be 



-2)/', 



and the temporal variation 



2'ipUc, 



^- "em 



a 
a 



(41) 



(42) 



If there were no scalar injection of energy and 
0, the Poynting theorem Eq. (|37p and the energy variation 
given by Eo. ip^ would lead to 



2lpUe 



2^e-^^ (B^+E^ ^_j.E + ^e-^^— (43) 
Stt 47r 



or 



j.E 



B2 

47r 



ipe' 



-2ip 



(44) 



As we consider phenomena where the motion of mat- 
ter is negligible, taking the first index as is equivalent 
to projecting along the fluid four-velocity. In addition, the 
total time derivative d/dt = d/dt + v.V will be equal to 
the partial time derivative d/dt. In the general case when 
there is viscosity and heat transfer, the right-hand side can 
be written, in the non-rclativistic limit, as 

- §iilpv' + ") + V.[pv(i^;2 ^ ^) _ ^ j]^ (45) 

where w is the specific enthalpy, u is the internal energy 
density, J is the heat flux, which can generally be writ- 
ten as — kVT. T is the temperature and k is the thermal 
conductivity. Finally, (v. a')k stands for Via'ii^, whe re a' is 
the viscous stress tensor, (jLandau fc Lifschitz Il987f ). As we 
stated above, we neglect the velocity of the fluid, so obtain 



t: 



(46) 



A note of caution regarding the internal energy is in 
order. We understand, as usual, "internal energy" as the 
energy that can be exchanged by the system in the pro- 
cesses considered (heat exchange, radiative transfer, etc.), 
which will differ from what wc understand by "rest mass" , 
which is the "non-convertible energy". If the scalar field 
can change the effective electric charge, then it can alter 
the electromagnetic contribution to the rest mass, and con- 
sequently, this contribution will be no longer "rest mass" , 
but "internal energy" . 

The time variation in the internal energy u will have 
two contributions: one corresponding to the cooling pro- 
cess ^Icooiing and another one related to the interaction 

with the scalar field This last term accounts for the 
dependence of the bulk of matter on the scalar field, which 
is mainly given by the electromagnetic contribution to the 
nuclear mass. Equation p2[) then will finally read 



du 



dt 



V- J 



B2 

An 



-2i/' 



47r 



(47) 



Since the scalar field is space independent, and given 
that the electromagnetic energy of matter is mainly ac- 
counted for the nuclear content, we assume that the follow- 
ing condition If - ^ 
obtain 



is fulfilled. We consequently 



du e-2V'B2?/; 

-T^UoUng + V • J « — . 



(48) 



This equation becomes clearer if we make a trivial 
change to produce 



V- J 



-^'AbV du. 



277 



dt 



I cooling ? 



(49) 



Kraiselburd et al.: Variation in a and energy production 



5 



which clearly shows that besides the standard cooling mech- 
anism of the body, there is a contribution from the partial 
release of the magnetic energy injected by the scalar field. 
Wc define 



a B2 



to be equal to twice the energy production per mass unit of 
any material substance a (using the approximation, e~'^'^ — > 
1 when << 1). 

We now consider our main physical assumption: the 
cooling term is not modified by the scalar field. The rea- 
sons for this assumption are fold: 1) as we have just shown, 
the electrostatic energy "injected" by the scalar field re- 
mains within the bulk matter (the cancellation of terms 
occurring as seen in Eg . (|48| ) . and 2) the thermal evolution 
should not change given the high thermal conductivity of 
the Earth and white dwarfs considered in this work. Thus, 
we expect the magnetic energy excess to be radiated away, 
increasing the heat flux J as shown in Eq. 



4. The electromagnetic energy of matter 

As we mentioned in the Sect.|31 the only "input" is that de- 
rived from the magnetic field. Stationary electric currents 
generated by charged particles and their static magnetic 
moments, and quantum fluctuations of the number den- 
sity are responsible for the generation of magnetic flclds in 
quantum mechani cs. These contrib u tions have been studied 
and calculated by iHauean fc Will ([l977h and MM (|l98l 
from a minimal nuclear s hell model using the following an al- 
ysis (for more details see iKraiselburd fc VucetichI (|2011f )'). 
The total magnetic energy of the nucleus can be written 

as 



^™ - ^ II / ^^^^ 



,{0\jix)\a)-{a\j{x')\0) 



\x — X' 



(51) 



where a runs over a complete set of eigenstates of the nu- 
clear Hamiltonian H. Neglecting the momentum depen- 
dence of the nuclear potential and assuming a constant 
density within the nucleus, we obtain the result 

{Q\j{x) \a) . {a\j{x') |0) ^ ^:^cos0, (52) 

where doa is the dipole density. Vat = is the nuclear 

volume, and 6 is the angle between x and x' . Hence, 



E,, 



i:aEL\do.\' Sdxdx'-^ 



2fi2c2 

f dxdx' , '^"^ ^, , 



y2 



(53) 



In the last equation, is equal to t-h— , and the 

first term can be computed from the connection between 
the strength function and the photoabsorption cross-section 



O'Oa = -^Eao\dao\ . 

he 

From this, we easily find that 

^ hc SEa{E)dE /• 

^E^,\d,^,\ ^-^-^-^ . j o{E)dE 



(54) 



(55) 



E / a{E)dE, 



where E ^ 25 MeV is the mean absorption energy, which is 
roughly independent of A (number of nucleons). 

The cross-section satisfies the Thomas-Reiche-Kuhn 
sum rule 

a{E)dE = (l-t-x) ^^ ^ ~ (H-a;)15MeVmbarnyl, 

mc A 

(56) 

where x ~ 0.2 takes into account exchange and velociy de- 
pendence of nuclear interactions. Combining Eas.([53|. ((55)) . 
and (|56p . we obtain 



E„.^ld^.^^±, /d3.rf3,.i(-)-^'(-') 



x — X' 



3 E 
20^ R{A)hc 



(57) 



adE, 



where R{A) is the nuclear radius. These quantities have the 
following approximate representation 

R{A) = 1.2Ai fm, J crdE ~ 1.6^ MeV fm^. (58) 

The fractional contribution of the magnetic energy to rest 
mass energy is then 



C(A) ^ 



.60465 X IQ-^A-^^^. (59) 



Table [T] shows typical values of C(^) computed using 
the semi-empirical mass formula and the contribution of 
neutrons and protons. 



Nucleus 




He^ 
Of 

Si?! 


5.42 
3.76 
3.41 
2.83 
2.25 



Table 1. ^ values for typical indoor stellar and planetary 
elements 



5. The Earth heat flux 

There are several models that attempt to explain the aver- 
age rate of secular cooling of the Earth in terms of varia- 
tions in the composition of the mantle melts through time 
(jLabrosse fc JaupartI l2007h . Constraints on these theories 
are set by terrestrial heat flow measurements on the sur- 
face. 

The contribution oi a/a to the heat flux can be calcu- 
lated using the global heat balance for the Earth, assuming 
that the machian contribution He is the only extra energy 
production, 



MeC\ 



dt 



— ^^tot 



+ Hc + Hc 



(60) 



where M e is the Earth's mass, Cp is the average heat capac- 
ity of the planet, is the mantle potential temperature, 
and Hq represents the heat generated by radioactive iso- 
topes. The total heat loss Qtot can be written as the sum 
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of two terms, one that comes from the loss of heat in the 
oceans Qoc, and the other fro m continental heat los s Qcon t- 
Using the results obtained by iLabrosse fc JaupartI (j2007l ) , 
we rewrite the total heat loss as Qtot ~ MCpXoTm, where 
Ag is the timescale constant for the secular Earth's cooling. 
Assuming that the most abundant elements of the Earth are 
oxygen, silica, and iron, the "extra" energy contribution can 
be written as 

He = Cc^Ho^, (61) 
where C is the mass-weighted averaged of the parameter 



Parameter 


Value 


Ho 


2.5 X 10"''' s"' 


Me 


5.94 X 10^* Kg 


Cp 


1200 J/Kg - K 


C 


2.75 X 10"'' 


Ag 


0.1 Gyr"^ 



Table 2. Values of parameters 



From Eq. (|14[) , we can describe the extra contribution as 
a function of time, writing ^lii as a power series, (jWeinberd 
[1973) 



a(t) 
ao 



1 + Hodt 



2 



(62) 



and then making a Taylor series expansion up to third or- 
der of He- Replacing this machian contribution in Ea. (|60p . 
solving it, and using the parameter's values from Table [51 
we find that the cosmological perturbation of the mantle's 
temperature AT™ in terms of the time interval At and 
is given by 



AT„(t) = 2.43 X lO^K/Gy, 
- 3.78 X 10^ K/Gy, 
-I- 3.05 X 10^ K/Gy 



a 



a 

a 
Hoa 



Hoa 



(At)= 



(63) 



At. 



According to ILabrosse fc Jaupar^ ()2007D . the total 
amount of cooling experienced by the Earth after an ini- 
tial magma ocean phase cannot exceed 200 K. Hence, in the 
past 2.5 Gyr, ATm < 200 K. With these restrictions, we 
obtain a bound for the time variation in a of. 



a 



< 1.93 X 10" 



Hoa 

Using this result in Ea.([T5|). we find that. 



< 0.15, 



< 0.39. 



(64) 



(65) 



A different bound can be obtained by observing that 
the total radiated power of the Earth Qtot can be ex- 
plained by radiactiv e deca y to within twenty per cent 
(jLabrosse fc Jauoartj (j2007D ). The most recent data was 



estimated from an adjustment made with 38347 measure- 
ments. The methodology was to use a half-space cooling ap- 
proximation for hydrothermal circulation in young oceanic 
crust; and in the remainder of the Earth's surface, the aver- 
age heat flow of various geological domains was estimated 
as defined by global digital maps of geology, and then a 
global estimate was made by r nultiplying t he total global 
area of the geological domain, (|Davies fc Davie s 20 nj). 

The result shows that Qtot « 47 TW ( in Table [S| this 
estimate is separated into continental and oceanic contri- 
butions). Therefore, 



\Qmach\ = \MECp\GT^{t)\ < 0.2Qtof 

In an interval of 2.5 Gyr, then we find that 



Hna 



< 3.98 X 10" 



(66) 



(67) 



and 



I <0.1, g<0.». 



6. The heat fluxes of the outer planets 

Jupiter, Saturn, Uranus, and Neptune are often called gas 
giants. They are massive planets with a thick atmosphere 
and a solid core. Jupiter and Saturn are composed primar- 
ily of hydrogen and helium, while Uranus and Neptune are 
sometimes called ice giants, as they are mostly composed of 
water, ammonia, and methane ices. By comparising the ob- 
served bolometric temperatures of giant planets with those 
expected when the planets are in thermal equilibrium with 
incident solar radiation, it is clear that all of these planets 
except for U ranus have a significant internal heat source, 
(llrwini 1200611 . 

In the case of Jupiter, the residual primordial heat emit- 
ted is caused by the continued cooling and shrinking of the 
planet via the Kelvin- Helmholtz mechanism. 

Saturn must have also started out hot inside like Jupiter 
as the result of its similar formation. But being somewhat 
smaller and less massive, Saturn was not as hot in the be- 
ginning of its life and has had time to cool. As a result, this 
planet has lost most of its primordial heat and there must 
be another source of most of its internal heat. This excess 
heat is generated by the precipitation of helium into its 
metallic hydrogen core. The heavier helium separates from 
the lighter hydrogen and drops toward the center. Small 
helium droplets form where it is cool enough, precipitate 
or rain down, and then dissolve at hotter deeper levels. As 
the helium at a higher level drizzles down through the sur- 
rounding hydrogen, the helium converts some of its energy 
into heat, (Lang 2003). 

The low value of Uranus' internal heat is still poorly 
undertood. One suggestion is that chemical composition 
gradients may act as inhibitors of heat transport from the 
planet's hot interior to the surface. Another hypothesis is 
that it was hit by a supermassive impactor that caused it to 
expel most of its primordial heat, leaving it with a depleted 
core temperature. Uranus has as much as 4Af^ of rocky 
materials hence, part of the internal flux (« 0.02VKm"^) 
comes from radioactive decay; Kelvin- Helmholtz mechanism 
would also be expected. 

Although Neptune is much farther from the Sun than 
Uranus, its thermal emission is almost equivalent. Several 
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Part of the Earth 


Area (lO'^m^) 


Heat flow (TW) 


Mean heat flow (21^) 


Continent 


2.073 


14.7 


70.9 


Ocean 


3.028 


31.9 


105.4 


Global Total 


5.101 


46.7 


91.6 



Table 3. Summary of the heat flow from lDavies fc Davied (|201Clf ) preferred estimates 



possible explanations have been suggested, including radio- 
genic heating from the planet's core, conversion of methane 
under high pressure into hydrogen, diamond and longer hy- 
drocarbons (the hydrogen and diamond would then rise and 
sink, respectively, releasing gravitational potential energy), 
and convection in the lower atmospher e that causes gravity 
waves to break abov e the tropopause (|Fortnev fc Hubbard 
[200llHubba7dl[l97l . 

We calculate the heat fluxes J^. for each planet using 
the equation of heat conduction 



dr \ dr J 



(69) 



where K is the effective thermal conductivity of the planet 
material. The heat flux is 



dr 



(70) 



If e is the planet mean heat production, which is esti- 
mated from the results in Table [1] according to the chemical 
composition of each planet, then 



dr 



'£(r')p(r')dr'=£-^, (71) 



7. Conclusions 



The energy exchange with ordinary matter in alternative 
theories with new fields such as Bcckenstein's theory is a 
quite undeveloped field of subject. Using the field equations 
and general hypothesis of the theory we have derived the 
energy transfer between matter and fields. The Hypothesis 
[5] in Sect. [5] is a probable key, because it states that the 
matter energy momentum tensor is the quantity that has 
to be added to the field sector in order to make the total 
tensor divergence free. We have also assumed that dark 
matter is electrically neutral, have neglected the motion of 
matter in the bodies considered, and have found that the 
dynamical feature of the electric charge makes the atomic 
electromagnetic energy part of the internal energy of the 
system. Eq. pS)) shows that there is another contribution 
to the heat current in addition to the cooling of matter, 
which is given by the time variation in the scalar field and 
the magnetic content of matter. We have also justified our 
assumption that the matter cooling rate is not modified 
by the scalar field. Finally, using a minimal nuclear shell 
model we estimated the magnetic energy content of matter, 
thus permitting us to quantify the anomalous heat fiux in 
terms of the fundamental parameters of the theory and the 
chemical composition of the body. 



hence the surface flux is 



3c = - K 



.dT_ 
dr 



= e 



m{R,) 



1 1 



(72) 



which is the fundamental equation. Thus, we compare 
the results of . with the observed fluxes obtained with 
Voyager 1, 2, and Cassini, (jPearl et al.lll990D . 

The "3(t" upper bounds and the corresponding 
"(^s/^p)" bounds are 



Planet 




6t 






Jupiter 


2.04 X 10"'' 


1.25 


Saturn 


1.55 X 10"^ 


1.09 


Uranus 


6.75 X 10"*^ 


0.72 


Neptune 


7.85 X 10"^ 


0.78 


Earth(^) 


1.93 X 10"^ 


0.39 


Earth'^' 


3.9 


? X 10"^ 


0.55 



Table 5. Bounds from the outer planets and the Earth C^^^ 
results from Egs. ljM)) and and from Eas. (|57|) and 

m)- 



Our tightest constraint was obtained by analyzing the 
geothermal aspects of the Earth, which are naturally the 
most clearly understood and reliably measured of our so- 
lar system, and the surface heat flux is very low. Our 
bounds are comparable with that obtained in the labo- 
rat ory by combining measure ments of the freque nces of 
Sr (iBlatt et all (l200l V Hg- H (iFortier et al.l (l2006h ). Yh+ 
(iPeik et al.l (I2004D) a nd H (iFischer et al.l (|2004D f relative 
to Caesium (iLi et al.ll201Cll : IUzanll2010l) . and only one order 
of magnitude weaker than Oklo's, which is the most strin- 
gent constraint o n a's time variation up to date (|Uzanll2010l : 
iFuiii et ani200Cl[ )). and another found from measurements 
of the ratio of Al -I- to Hg-|- optical clock frequencies ov er a 
period of a year (jRosenband et al.ll2008t iLi et al.l[2010h fsee 
Table |6]). The constraints we found depend on the cool- 
ing model of the Earth, but the re is a genera l agreement 
about the mechanisms behind it (jjessopl (|l99Cl[ )l The outer 
planets have provided us with additional constraints, which 
are between the same and one order of magnitude weaker 
than Earth's, but are nevertheless valuable, as the chemical 
composition and cooling mechanisms differ widely from our 
planet. The data set considered here is able to place inde- 
pendent constraints on the theory parameters. This analysis 
may be applied to other theories with extra fields that in- 
troduce extra "internal energies" to matter. We will report 
on further work in future publications. 
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Kraiselburd et al.: Variation in a and energy production 



Planet 


Jobs (W/m^) 


M (Kg) 


R (m) 


J<. (W/m^) ^ 


Jupiter 


5.44 ±0.43 


1.90 X 10^^ 


7.14 X lO'^ 


6.35 X 10-* 


Saturn 


2.01 ±0.14 


5.68 X 10^*^ 


6.03 X lO'^ 


2.71 X 10^ 


Uranus 


0.042 ± 0.047 


8.68 X 10^^ 


2.556 X lO'^ 


2.08 X 10^ 


Neptune 


0.43 ± 0.09 


1.02 X 10^** 


2.47 X lO'^ 


3.44 X 10^ 



Table 4. The observed heat flux, mass, radius, and calculated heat flux of the outer planets 



Constraint 


If 1 (yr-^) 


Reference 


ClocksCs 


(3.3 ±3.0) X 10""" 


(1) 


ClocksHg 


(5.3 ±7.9) X 10"^'^ 


(1) 


Oklo 


(2.50 ±0.83) X 10""^ 


(2) 


Je 


1.52 X lO-""^ 


(3) 




3.14 X 10""'^ 


(4) 


J J up 


1.61 X 10"^^ 


(5) 


J Sat 


1.22 X 10""^ 


(5) 


Jur 


5.32 X 10"" 


(5) 


Jncp 


6.19 X 10"" 


(5) 



Table 6. The table comparises different kinds of con- 
straints, the value of — , and the reference. References (1) 
iLi et all (|20T1 : (2 lFuiii et al.l (l2000l) : (3) EqlMl (4) EqEH 
(5) Tahl^ 



Uzan, J. P. 2003, Rev. Mod. Phys., 75, 403, hep-ph/0205340 
Uzan, J. P. 2010, arXiv:1009.5514vl [astro-ph.CO] 
Weinberg, S. 1972, Gravitation and cosmology (New York-USA: John 
Wiley and Sons) 

Will, C. M. 1981, Theory and Experiment in Gravitational Physics 
(Cambridge- London: C. U. P.) 
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